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Abstract 

It is demonstrated that in simple soliton models essential features of the 
electro-magnetic nucleon form factors observed over three orders of magni- 
tude in momentum transfer t are naturally reproduced. The analysis shows 
that three basic ingredients are required: an extended object, partial coupling 
to vector mesons, and relativistic recoil corrections. We use for the extended 
object the standard skyrmion, one vector meson propagator for both isospin 
channels, and the relativistic boost to the Breit frame. Continuation to time- 
like t leads to quite stable results for the spectral functions in the regime from 
the 2- or 3-pion threshold to about two rho masses. Especially the onset of 
the continuous part of the spectral functions at threshold can be reliably de- 
termined and there are strong analogies to the results imposed on dispersion 
theoretic approaches by the unitarity constraint. 
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I. INTRODUCTION 



Topological soliton models for structure and dynamics of baryons are based on effective 
nonlinear lagrangians for selected mesonic degrees of freedom. These usually comprise the 
pseudoscalar Goldstone boson octet of spontaneously broken chiral symmetry, but also the 
light vector and axial vector mesons have been included. 

A decisive advantage of the soliton concept as compared to models where explicit point- 
like fermion fields are coupled to meson and gauge fields is the fact that already in leading 
classical approximation the spatial structure of the baryon as an extended object is obtained 
from the underlying effective action. Therefore all types of form factors can readily be ex- 
tracted from the models, and comparison with the experimentally observed dependence on 
momentum transfer presents a stringent test for the resulting spatial profiles. Specifically, 
the electro-magnetic form factors of the nucleon for which we expect a wealth of precise data 
in the few (GeV/c) region from the new generation of electron accelerators pose a severe 
challenge for chiral soliton models. 

Shortly after the initial work of || the e.m. nucleon form factors have been evaluated for 
various versions of effective meson models and various sets of parameters, mainly for 
momentum transfers Q < 1 GeV/c. The conclusion was that the resulting Q 2 - dependence 
follows roughly the standard dipole form, and that reasonable values for radii require the 
presence of explicit vector mesons. Unfortunately, in cases where these results have been 
included in more detailed comparisons with new data (e.g. in P|,|TU||), they have left the 
prejudice that soliton form factors are not satisfactory. However, for momentum transfers 
Q > 1 GeV/c it is important to incorporate relativistic kinematical corrections into the 
form factors. This is crucial for quark bag or cluster models |)-10] as well as for soliton 



models ||11|| . The implementation of these corrections is fairly easy for solitonic nucleons due 
to the Lorentz covariance of the underlying field equations. 

The proton magnetic form factor is quite accurately known up to 30 (GeV/c) 2 and it is 
not difficult to show [F2 that with suitable choice of coupling constants the relativistically 
corrected soliton form factor calculated from a minimal vector meson model follows the data 
closely over more than three orders of magnitude in Q 2 . This may be a bit disappointing 
because one might have hoped that in this region of high momentum transfer the character- 
istic shape of form factors could reflect dynamical features of QCD 113]. On the other hand, 
a similar conclusion had already been reached by Holder |14| through a dispersion-theoretic 
fit of the spectral functions which included a continuous part to account for the Frazer-Fulco 



unitarity constraint ||T5 |. 



In the following we would like to analyze which basic qualitative features of the e.m. 
form factors and their spectral functions can be expected from the soliton concept for nu- 
cleons and to which extend these may correspond to observed facts. In order to expose the 
typical features of solitonic e.m. form factors in the most transparent way it is necessary to 
choose the essential ingredients as simple as possible. We therefore do not strive for optimal 
agreement with the data for some sophisticated model, but instead choose the most simple 
construction: i) a purely pionic soliton created through the standard Skyrme term; this 
represents the extended object with its spatial structure defined through the chiral profile 
.F(r); ii) its minimal coupling to the electro-magnetic field is partly mediated through vec- 
tor mesons of mass m v represented by one Klein-Gordon propagator; in) relativistic boost 
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factors to the Breit frame account for kinematic corrections. 

In section II we show that these three ingredients are sufficient to produce for Q 2 < 1 
(GeV/c) 2 perfect agreement with the standard dipole and for Q 2 > 1 (GeV/c) 2 the charac- 
teristic deviations from the dipole form. We also get indications at which point more precise 
data could lead to discrepancies. 

Of course, it is most interesting to try to learn something about the spectral functions 
which underlie the resulting form factors. Especially intriguing is the question if and what 
the Skyrme soliton knows about the Frazer-Fulco constraint. A direct investigation of this 
problem by relating the pion-baryon phase shifts to the timelike form factor appears difficult, 
because the 7r-iV S-matrix in the soliton formulation depends only on one variable. Earlier 
analyses of the form factors for timelike momentum transfers fU are not helpful in this re- 
spect because the method of analytic continuation used there fTEfl , provides only averaged 



form factors with an averaging interval that conceals all detailed structures. Unfortunately, 
even within this simple model the continuation into the timelike regime remains ill-posed; 
but it seems that obtaining the spectral functions as Laplace inverse of the spatial densities 
provides rather stable results in the most interesting region from threshold to about 2 m g . 
We discuss this procedure in detail in section III for the electric form factors; the magnetic 
spectral functions then are easily obtained by differentiation or integration. Once the spec- 
tral functions of the plain skyrmion are determined in its restframe, the modifications due 
to partial coupling to vector mesons, and the boost to the Breit frame pose no additional 
problem. They are discussed in sections IV and V. 

II. RELATIVISTIC FORM FACTORS IN THE SOLITON MODEL 



We create the nucleon as topological soliton through the standard Skyrme model (T? 
(U e SU(2)) 

C(U) = C {2) + £ (4) (2.1) 

f 2 r 

£( 2 ) __ ±ii_ I (-trL^ + mltr(U + U ] - 2)) d 3 x. (2.2) 

C^ = ^Jtr[L„L v ] 2 d 3 x (2.3) 

where denotes the chiral gradients 

Lfj, = U^d^U. (2.4) 

The pion decay constant f n and the pion mass m n have their physical values f n =93 MeV 
and m 7r =138 MeV, the Skyrme parameter e is taken at e=4.25 which leads to the correct 
Delta-nucleon mass difference of Ma — M^=29b MeV. 

The static soliton solution, the Skyrme hedgehog U = exp(if-fF(r)) is characterized by 
the numerically determined chiral profile F(r) which carries all the information about the 
spatial structure of the extended object. 
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In the Breit frame for spacelike t = —Q 2 the isoscalar (T=0) and isovector (T=l) electric 
and magnetic (Sachs) form factors are given as Fourier transforms of baryon density B Q (r) 
(for T=0) and moment-of-inertia density Bi(r) (for T=l): 

G T E =°(-Q 2 ) = \j d*r Jo (Qr)B (r) (2.5) 
Glf{-Q 2 ) = Ald?r # r 2 B (r) (2.6) 



2r 2 B J Qr 

G T E =1 (-Q 2 ) = \J d 3 r Jo (Qr)B 1 (r) (2.7) 
Glr\-Q 2 ) = lJd\ 3 ^B 1 (r) (2.8) 



Expressed in terms of the chiral profile F(r) the densities Bo(r) and Bi(r) are 

1 pi ®n 
2tc 2 r 



^r) = ^F'^, (2.9) 



^ = Se- 



isin 2 F+4sin 2 F( F' 2 + 



sin 



2 F 



(2.10) 



with normalization / d 3 r Bq(t) =1, and / d 3 r B\(r) =1. 

For the magnetic form factors the isoscalar and isovector magnetic moments 



2/io = ftp + A*r. 



M 



2 



30 



2/ii = n P - jj n = -M N Q 



(2.11) 



(with baryonic square radius r 2 B = J d 3 r r 2 B (r)) have, respectively, been divided out such 
that they both satisfy the same normalization condition as their electric counterparts 



G° 



M 



(0) 



g\M = -. 



(2.12) 



The Skyrme term in (|2.3|) effectively represents the lowest local approximation to 
explicit inclusion of higher resonances, notably g and uj mesons [|T8|-f20|j. Propagating vector 



mesons introduce additional pole structures near the vector meson mass into the form fac- 
tors, which cannot be accounted for by the local Skyrme term alone in a satisfactory way. 
However, instead of dealing with the ambiguities and numerous parameters of vector meson 
models it appears more appropriate to incorporate their effect into one common factor A(t) 
to be multiplied with the pure Skyrme model form factors 



D-(P 



A(t) = A 



mi 



+ (1-A). 



(2.13) 



With the physical vector meson mass my=770 MeV this factor contains one additional 
parameter A which allows to interpolate between complete vector meson dominance (A=l) 
and the purely pionic Skyrme model (A=0). The resulting form factors 
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G°E 1 , M (t)=A(t)G^f I \t) (2.14) 

can be shown to represent good approximations to results obtained in models with explicit 
vector mesons included. (For the coupling to the a;-meson in the minimal model it is exact.) 



Although the definitions (|2.5|) - (|2.8|) hold in the Breit frame only, expressions (|2.9| ) and 
( p,10| ) for the densities contain the chiral profile F(r) as obtained in the soliton restframe. 
However, both frames coincide only for infinite soliton mass M$. For values of Q 2 near and 
beyond the actual soliton mass Ms kinematic corrections have to be included by employing 
the hedgehog soliton boosted to the Breit frame, where it moves with velocity v 

7 2 = (1 - vY 1 = l + -^- = l--L- (2.15) 

1 V ; 4M 2 AM 2 V ; 

with spacelike momentum transfer t = —Q 2 < 0. 

An approximate way to implement these recoil corrections has been suggested by Ji |TT . 
The electric form factors are Lorentz scalars, therefore one has 

G E {t) = Gt\th 2 ) (2.16) 



where is the non-relativistic form factor evaluated in the soliton restframe. 

For the magnetic form factors the boost to the Breit frame introduces an additional 
factor 7 -2 

Guit) = - 2 G&\t/y 2 ) . (2.17) 

7 

These boost factors are well-known in models for relativistic quark clusters |J. There, 
however, Ge carries the same factor 7~ 2 as Gu due to the normalization of the spectator 
quarks. 

The kinematical transformations fl2.16|) and (|2 . 1 7|) determine the asymptotic behavior of 
the relativistic form factors for large spacelike t — > — oo: 



G E (t) 



A A/f 2 fl 

G% r \-AM 2 ) + ^^G^ r \-4M 2 ) 



+ 0[±) (2.18) 



G M (t) 



1 



2 U M I 



-4M 2 



+ o 



r 2 



(2.19) 



These asymptotic forms show a rather undesirable feature of the boost transformations 
Q2.16j ), Q2.17D because they violate the super convergence rule expected for e.m. form fac- 
tors IB1IE2I 



lim tG EM (t) = 0. 
t— >— 00 



(2.20) 



Practically, however, the absolute values of G^ r ^ and its derivative at t = —AM 2 are so 
small that up to Q 2 =10 (GeV/c) 2 the influence of the undesired terms is not important. 
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This is demonstrated in figs. 1,2 where the dashed lines show the e.m. form factors after 
subtraction of the terms written out in square brackets in eqs. (|2.18| ), ( [2.19| ). 



In principle, within tree approximation, the kinematic mass M must be identified with 
the classical soliton mass Mg. Ideally, of course, M should coincide with the physical 
nucleon mass M/v. However, it is known, that this difference is related to quantum cor- 
rections [^,0] which probably also affect the existence and positions of zeros in the non 



relativistic form factors. Therefore, M is not really well defined in the model and could be 
used as additional parameter to minimize the undesired terms in ( |2.18|) , (|2.19|) by choosing 
M such that G { E r) M {-AM 2 ) = 0. In this way superconvergence could be enforced at least for 
Gm- However, to keep the number of parameters as small as possible we strictly adhere to 
the tree approximation and identify in the following the kinematical mass M in the boost 
transformation with the soliton mass M$ which is Ms=1648 MeV for the Skyrme model 
with e=4.25. 

Thus the only parameter remaining in the model is A. It enters sensitively into the radii 
and interpolates between the purely pionic radii (which are generally too small) and the 
completely vector meson dominated radii (which are too big): 

<r2> =7^ 7T+ G W = X A + < r * > Pl ° mC ■ ( 2 - 21 ) 
G(0) at m v 

(For the electric neutron radius we take 

< r 2 > E = 6 ^G E (0) = ~(< r 2 >%-< r 2 >%) ; (2.22) 

here the vector meson contributions cancel). 

The form factors shown in figs. 1,2 are calculated with A=0.75, i.e. the e.m. coupling 
is strongly but not completely vector meson dominated. The resulting square radii are (in 
fm 2 ): 

r\ = 0.233; < r 2 > P E = 0.756; < r 2 > P M = 0.717; < r 2 > n E = -0.228; < r 2 > n M = 0.746; 

Up to Q 2 ~ 0.5 (GeV/c) 2 , G P M , G P E , G n M follow the standard dipole G D (t) = (l-t/0.71)- 2 
very closely. The magnetic form factors rise above the dipole in the region between 0.5 and 
5 (GeV/c) 2 , and fall below the dipole above 5 (GeV/c) 2 . The neutron form factor G\ { rises 
slightly above the proton G P M , which seems to be in conflict with the data of |^5| for G 1 ^. In 



fig. 2 for Gyi in addition to more recent data points [^5H 28|| we have also included a collection 



of older data [p9j-p4| . They seem to indicate significant deviations from the standard dipole 
parametrization also for small Q 2 . But we find that at least up to 1 (GeV/c) 2 the ratio 
G^/Gm stays constant with good accuracy, and it appears difficult to accommodate such 
differences between proton and neutron data. With explicit dynamical inclusion of vector 
mesons it is not difficult to find a parameter choice such that the rise of G P M above the dipole 
coincides closely with that observed in the data (see e.g. |I^|). but it should be remembered 
that both, isoscalar and isovector magnetic moments /i and /ii are known to be subject to 
sizable quantum corrections [23J], which may also affect the shape of the form factors. 

In G P E the corresponding rise above the dipole is suppressed and the form factor is 
dominated by a zero near 10 (GeV/c) 2 . From our discussion of the boost transformation 
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it is clear that the high-Q 2 behavior of the ratios G/Gd is sensitive to the position of the 
first zeros in the form factors (relative to AM 2 ). From ( |2.5|) - (|2.8|) we expect that the first 
zero in Gg occurs at lower Q 2 than the first zero in Gm- (It may be noted that in the 
plain Skyrme model Gm is monotonously decreasing and has no zeros; however, inclusion 
of sixth-order terms or explicit vector mesons leads to small fluctuations and corresponding 
zeros also in Gm)- As long as these conditions are not inverted by quantum corrections we 
would therefore conclude that the rapid fall-off in Ge relative to Gm is a typical and rather 
stable feature of soliton models. A similar behavior is, however, also observed in relativistic 
quark bag or cluster models (see e.g. f7|,PJiT| ) . 

Data points for the electric neutron form factor extracted from e-d scattering depend 
sensitively on the choice of the deuteron wave functions. This is illustrated in fig. 2 for G E 
by the dotted lines which represent fits of |35[ of the form 



G E (-Q 2 ) = -a^ n ^G D {-Q 2 ){l + b^)- 1 (2.23) 

where a and b depend on the choice of the n-p potential used in the deuteron wave function. 
Included in fig. 2 is the data set of |3f| for the Paris potential and an older data collection 
of IBBI based on the Feshbach-Lomon wave functions. The dash-dotted line is the Galster 



fit [36| with a—1, 6=5.6. The shape of the calculated G E (full line) follows closely the 
Galster parametrization. (The absolute values of G% and thus the value of the slope near 
Q 2 = can easily be reduced by allowing A to be slightly different for isoscalar and isovector 
mesons. Again, however, without quantum corrections it is not very meaningful to adjust 
parameters to the very sensitive electric neutron square radius, which experimentally is near 
< r 2 > E =-0.12 fm 2 ). 

So, apart from choosing A=0.75 we have made no attempt to further improve agreement 
with the experimental data. The results sufficiently demonstrate that this extremely simple 
model is capable to cover the essential features of the observed form factors. Significant 
quantum corrections are expected to the magnetic moments, to the radii, and to M. By the 
choice of A we have compensated their possible influence on the radii. The high-Q 2 behavior 
of the form factors is sensitive to the choice of M and expected differences in loop corrections 
could possibly be absorbed into slightly different values of M for different form factors. But 
here all form factors in figs.l and 2 have been calculated with the same value of M (the 
soliton mass). As long as rotational and loop corrections to the currents are not included, 
the accuracy shown in figs. 1,2 is satisfactory and we now turn to the typical features of the 
spectral functions which underlie these form factors. 



III. SPECTRAL FUNCTIONS FOR NUCLEON FORM FACTORS IN THE 
NONRELATIVISTIC SKYRME MODEL 

It is generally assumed that e.m. form factors G(t) fulfill unsubtracted dispersion rela- 
tions 

1 r°° r(t') 

G(t) = - -^-dt> (3.1) 

7T Jt t — t 



7 



where the spectral function F(t) = Im G(t) is given by the boundary value of G(t + iO) along 
the upper edge of the cut which extends from i along the positive real axis in the t-plane. 
The lower limit is to — 4m 2 for isovector and to — 9m 2 for isoscalar form factors. The de- 
termination of spectral functions from experimental form factor data in the spacelike region 
t < requires an analytical continuation to the cut which is known to be mathematically 
an ill-posed problem. Therefore such extrapolations have to be performed under certain 
stabilizing assumptions and the resulting spectral functions reflect these assumptions. 

It is therefore of interest to find out what successful nucleon models may tell us about the 
structure of the spectral functions. But it should be noted that (unless the models provide 
analytical expressions for the form factors) in practice we face the same ill-posed problem of 
inverting (|3.1| ) with G(t) given with finite numerical accuracy for all values t < to- However, 
apart from the fact that the numerical errors of the 'model data' for G(t) can be made 
(arbitrarily) small, the models in general will provide at least some analytical constraints on 
form factors (or on the underlying spatial densities) which will help to stabilize the analytic 
continuation. 

In the following we analyze this procedure for the Sachs form factors of the simple Skyrme 
model. Inclusion of additional explicit vector meson poles and relativistic corrections then 
will be straightforward. 



A. Isoscalar electric G° E {t) 

We repeat that in the Breit frame for spacelike t = —Q 2 the isoscalar electric form factor 
is given by the Fourier transform of the baryon density B (r) 

G° E (-Q 2 ) = \j d ' r UQr)B (r) . (3.2) 

If we assume the existence of a function T(t') such that the integral in ( |3.1| ) is convergent 
then spatial density B (r) and spectral function T(t') are related through 

rB (r) = — / e-^lV)^ (3.3) 

(with notation t' = /i 2 ), i.e. the function fiT(n 2 ) is the Laplace inverse of 7r 2 rB (r). Although 
mathematically the problem (|3.3|) is no less ill-posed than fl3.1|) , models generally supply 
analytical information for the behavior of Bo(r) for r — > and r — > oo, in which case ( |3.3| ) 
appears better suited for obtaining the spectral functions. 

Specifically, the baryon density Bq(t) for the Skyrme hedgehog is given by 



F' sin 2 F 



and the asymptotic form of the chiral angle F(r) is 



F(r) ™ ^r(l + m^r)e~ m - r (3.5) 
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and F{r) ^ n - Br + • • • (3.6) 

where the constants A, B are fixed by solving the nonlinear equation for the soliton profile. 
The strategy for constructing r(yU 2 ) then is to approximate the numerically determined 
density rB (r) by a set of functions with known Laplace inverse such that the constraints 
imposed by fl3.5|), ( |3.6| ) are fulfilled and ( |3.3| ) is satisfied for all values of r. 

Commonly, spectral functions are set up as sums of several discrete 'monopoles' at po- 
sitions fi = vq with strengths ao 

r^VH^ao^/i 2 -^ 2 ). (3.7) 

Similarly, one may consider sums of 'dipoles' at positions fi = v\ with strengths a\ 

r (1) (^ 2 ) = E«i^V-^i 2 )- (3.8) 

Their contributions to rB^ir) according to ( |3.3| ) are of the form 

2^rB (r) = a e'^ + J> (?f) e"^} . 
In addition to fl3.7|), (|3.8|) we define higher-order pole structures in the spectral function 



T (2 V) = a 2 (z, 2 VV - 4) - - vl)) , (3.9) 

r (3 V) = «3 (V'V - ^ 2 ) - - , (3-10) 

I^V) = a 4 (^"V - v\) - 3z/ 4 6 r(^ 2 - z/ 2 ) + - i/ 2 )) , etc. (3.11) 

such that their contributions to rB^ir) are 

2^rB (r) = { £ a,(^)Y-'}. (3.12) 
i=o,i,.. v A 7 

This convenient form proves very efficient for an accurate representation of numerically 
obtained functions rB (r); we therefore prefer to use the pole combinations rw instead of 
the 'pure multipoles' 

Evidently, however, the representation (|3.12|) does not allow to satisfy the asymptotic 
constraint ( p.5|) which requires that rB (r) decreases for large r like 



2n 2 rB (r) - A z ^-e~^ r (3.13) 

with degree-four polynomial 

N 4 (r) = (1 + m w r) 2 (2 + 2m w r + (m w r) 2 ) . (3.14) 

In order to correctly reproduce the asymptotic behavior of the spatial density a term like 
( |3.13j ) has to be included on the right-hand side of fl3.12| ). However, to be able to satisfy the 
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constraint (|3.6|) for r — > 0, the denominator in ( |3.13|) has to be replaced by a degree-eight 
polynomial 

8 

D s {r) = Y,c i r i with eg = 1. (3.15) 

i=0 

Then, the analytical knowledge about rBo(r) and its derivatives at r = fixes the amplitudes 
of monopoles, dipoles,... in terms of the coefficients Cj. In contrast to the discrete structures 
( p.7|) - (p.ll|) of the spectral function which underlie the right-hand side of (|3.12| ) the rational 
function N±(r) / D s (r) corresponds to a continuous part r cont (/i 2 ) of the spectral distribution 
r(/x 2 ) which (due to the exponential factor exp(— 3m n r)) is nonzero only for \i > 3^71^). Its 
form near threshold is determined by the asymptotics fl3.5| ). 

Through an extended analysis of the nonlinear differential equation for F(r) additional 
terms in (|3.5|) for r — > oo can be determined which lead to improved numerator polynomials 
Nk(r) in ( p,13| ) with k > 4. Further terms in the expansion of sin 2 F in ( |3.4| ) carry the 



exponential factor exp(— 5m n r) and create additional continuous contributions along the 
cut for > 5m n . The evaluation of further terms of F(r) in (|3.6|) for r — > provide 



additional equations to fix more amplitudes in ( 3.12|) . In this way increasingly accurate 
approximations to rB (r) could be constructed. 

There is, however, a severe constraint which has to be imposed on acceptable denominator 
polynomials D^(r): the Laplace inverse of Ni/D^ (where k > i) contains the exponentials 
exp(«j/x) where «j are the roots of Dk- Therefore the integral in (|3.3|) converges only for r > 
Re a max , where a max is the root with the largest positive real part. Consequently all roots 
of acceptable denominator polynomials must lie in the left complex half plane. (For simple 
roots the imaginary axis is still acceptable). 

The fact that the "data", the numerically obtained function B Q (r), is known with great 
precision allows to create very satisfactory fits by including a sufficient number of terms at 
different positions z/j in the sum on the r.h.s. of 



2^rB (r) = A^e~^ + { £ * )V*'} (3.16) 



'i=0,l, 

and constrain their amplitudes di through boundary conditions at r = 

rB (r)\ = 0, (rB (r))% = B (0) , (rB (r))"\ = , (rB (r))'"\ = 3£ '(0) , (3.17) 

or, explicitly, 



= 2A 3 — + a (3.18) 



Co 



2n 2 B (0) = -2A 3 % - a u + -a 1 v 1 
c l 



= -4A 3 m ' C ° + C 3 ° C2 ° { + a,vl - *A + \*A 
eg 2 



6^(0) = 2A^ A + 6m ' C ° Cl + U A C0ClC2 ~ 6C ° C3 ~ 6C? - a ul + 2 fll *? - ^ + 3 3 

Cn 2 2 4 
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The first of these equations shows that at least one discrete monopole is necessary to allow 
for rB (r) | = 0, and that its discrete monopole strength a , however, is completely com- 
pensated by the continuous part of the spectral function, because we have from ( |3.16| ) and 
XI that 

TcontGuW = A 3 ^\ r=0 = —. (3.19) 

So, evidently, the origin of the dipole nature of the form factor lies in the fact that the 
density Bo(r) has a finite limit at r=0. 

Due to the fact that the set of functions which appear on the r.h.s. of ( |3.16| ) is neither 



orthogonal nor complete, there remain severe ambiguities in their selection. From the phys- 
ical point of view we might prefer to include only a set of monopoles at different positions. 
This, however, proves very inefficient for an accurate fit. Including more terms in (|3.16|) 



leads to improved fits; however, as usual, it does not make sense to include many more 
terms, because different fits of comparable (improved) accuracy may differ appreciably in 
their resulting distribution of multipole strengths. We find the best results (with respect to 
stability and accuracy) by including the four terms with i = 0, ..,3 in ( |3.16| ), i.e. just one 
term for each value of i at positions v\ with the four amplitudes <2j fixed through the four 
relations ( p. 18 ). The resulting pole structure can only be interpreted as an effective discrete 



representation of some undetermined underlying strength distribution. As is evident from 
( p,18| ), location and strength of these effective pole structures are in close connection to the 
additional continuous part which is severely constrained through the asymptotic behavior of 
the spatial density. Thus (|3.16| ) with i = 0, .., 3 appears as a simple and natural form which 
satisfies the asymptotics for r — > and r — > oo and allows for a very accurate fit of rBo(r). 

Allowing the roots of D$ to vary freely in the left complex half plane leads to two distinct 
double roots and their complex conjugates, i.e. 

D 8 = \(r-a 1 )\r-a 2 ) 2 \ 2 . (3.20) 

Imposing this form for D$ leaves two complex roots ai, a 2 and four pole positions z/j, 
{i = 0, .., 3) as parameters for the fit (|3.16|) . We find that the position of the roots ai, a 2 is 
quite stable with respect to different pole combinations, and the smallest \ 2 is obtained for 
the combination given in ( |3.16| ). 

With (|3.20| ) the decomposition of N^/Dg into partial fractions is 



Ni{r) 7i 7u 72 722 . foot\ 

- T7 — = + - + + - ^ + cc. (3.21) 

L> 8 {r) r — ai (r — a\) z r — a 2 (r — a 2 ) z 

where the coefficients 7j satisfy relations like 

Re ]T 7i = 

i=l,2 

Re (7« + a*7i) = 

i=l,2 

Re Yl (2«i7ii + ahi) = (3.22) 

i=l,2 
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and so forth, which imply that the Laplace inverse 7(/i) of N^/Dg 



7^) = 2Re (7, + T« (3.23) 

t=l,2 

satisfies at threshold \i = 

7 (0) = 0; 7 '(0)=0; T "(0) = 0; t"'(0) = mj; etc. (3.24) 

(These relations, of course, follow also quite generally by comparing the coefficients for 
positive powers of r for r — > oo in the definition of 7(/i) 



Ni(r) = D k (r) e'^j^dfi (3.25) 



for polynomials iVj(r) and Df.(r).) From ( |3.3|) and (|3.16|) the continuous part of r(/x 2 ) then 
is given by 

A 3 

r c0 nt(/U 2 ) = 7^70" - 3m,,) (3.26) 
for /i > 3m^. At threshold fi 2 — > (3m ff ) 2 its behavior is fixed by (|3.22|) as 

r c ont(/i 2 ) -> - (3m^) 2 ) 3 . (3.27) 

It should be noted that this onset of the spectral density at threshold is completely de- 
termined by the form of the baryon density (|3.4j ) and the asymptotics of the chiral angle 
( P-5| ) and is independent of the discrete structures added in (|3.16| ). It cannot be obtained 
in parametrizations of spectral functions in terms of discrete pole structures only, but the 
existence of this continuous part affects positions and amplitudes of poles through equations 
like (|3T8D . 

Altogether, the isoscalar electric spectral function is obtained in the form 

r /=0 (t) = r cont (t) + ^r«(t) (3.28) 

where the continuous part is given by ( |3.26| ) with ( |3.23j ) and the amplitudes of the discrete 
pole structures (|3.7Q - (|3.10|) are fixed through (|3.18|) . 

For jj? ^> (3m n ) 2 the continuous part r cont (/i 2 ) is dominated by the exponential expaj/i 
corresponding to the root a, = — £j ± itpi of D$ with the smallest (absolute) value of its real 
part £j, i.e. r cont oscillates with slowly decreasing amplitude. It turns out that the imaginary 
part fi (i.e. the period of the oscillations) is sharply fixed by the fit, while variations of the 
real part — within a limited range affect the quality of the fit only very little. (By putting 
a constraint on £j the oscillations thus can be damped without severe consequences for the 
fit). 

As an example we present the result of such an analysis for the standard Skyrme model 
(with f n = 93 MeV and e = 4.25). A typical fit leads to the (two-fold) roots (in units of 
(inverse) m g = 770 MeV) 
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ai = -0.725 ±i 1.505 ; a 2 = -0.5 ± % 1.823 [m" 1 ] 



(3.29) 



Here the real part of «2 has been kept fixed at -0.5. If e 2 is allowed to move freely it will 
slowly approach ^2=0 with only marginal improvement of the fit and little change in pole 
positions. The choice ( |3.29| ) efficiently cuts down the oscillations in the continuous part of 
the spectral function for large \x (see fig. 3). The pole locations and strengths for the same 
fit are listed in table I. 





i = 


i = 1 


i = 2 


i = 3 


Ui [m e ] 


1.61 


5.14 


1.93 


2.09 


Oj [m 2 Q ] 


-8.80 


0.47 


15.99 


-4.89 



TABLE I. Positions and strengths of the discrete structures for the isoscalar electric spectral 
function (in the soliton restframe) 



A very stable feature of this result is the discrete monopole strength T^(t) near 1.6 m g . 
Of course, this discrete monopole strength is completely compensated by the continuous 
part r cont (t) which is very small near t « m 2 Q but has a first pronounced maximum near 
(2.5m g ) 2 (see fig. 3). Also this feature of r cont (t) is very stable. There is no discrete low- 
lying dipole strength. But contains dipole strength and strong quadrupole strength 
near 1.5 GeV. Together with some octupole strength this may indicate broader underlying 
structures, although the details are certainly peculiar to our specific choice of allowed pole 
structures in the ansatz ( |3.16| ). Altogether it appears that in this way we obtain a quite 
reliable picture of the spectral function which underlies the isoscalar electric form factor of 
the Skyrme model up to about 2-3 m g . 



B. Isovector electric G\(t) 
For the isovector electric form factor 

GU-Q 2 ) = \j d " r MQr)Bi(r) (3.30) 
the baryon density B (r) in (|3.2|) and ( |3.3|) is replaced by the (isorotation) inertia density 

Bl (r) = ^F(H + + ^)) , (3.31) 

normalized by the moment of inertia G such that / Bi(r)d 3 r = 1. The asymptotic constraint 
( |3.5|) then suggests a rational function N2/D4 on the right-hand side in the representation 
of Bi{r) (corresponding to (|3.16|) ) 

J^BiCr) = h^f\e~^ + -{E«oe-^ + ^^'^ + ' ' ) ( 3 - 32 ) 

with 
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4tt 2 

h = ^f 2 A 2 

N 2 (r) = (1 + m^rf 

4 

D A = Y^c i r i with c 4 = 1 . (3.33) 

i=0 



The boundary conditions at r = are 



4 _o / „o 2 



B x (0)=0; B[(0) = 0; B'l(0) = ^B 2 (f 2 + -B 2 ); B»'(0) = 0; (3.34) 

These conditions exclude discrete monopole strength. It turns out that at least four discrete 
structures are needed for satisfactory fits. Again we find that by far the best result is 
obtained if we allow for just one term of each power of r in ( |3.32| ), i.e. we use 

2n 2 Bl (r) = h^-e-^ + ^E^y)' 6 ""} > t 3 ' 35 ) 
where the four amplitudes dj are fixed through relations ( [3.34 ), which read explicitly 

c 2 

„ Ci v 2 v\ 

= _ 2/l M^iijp^i) + a »l _ a2 | + „ 3 i 

, Amlcl + 6m^cicg - 6cgc 3 + 12c CiC 2 - 6c? v\ 3z/| 3z/^ 

= U — - — j-^ — ai— + a 2 — - — a? — - + a 4 — - 

J 4 2 4 3 4 8 

The two complex roots aci, a 2 of the denominator polynomical -D 4 (r) of degree 4 

D 4 (r) = |(r- ai )(r-a 2 )| 2 (3.36) 

together with the pole positions Uj (z=l,..,4) provide 8 real parameters for a very accurate 
fit. Again the roots a« = — £j ± are restricted to the left complex half plane. 

It is instructive to notice that (in contrast to the isoscalar case) the different asymptotic 
behavior of Bi(r) for r — > oo causes a discontinuity of the spectral function at threshold 
fi 2 — > (2m 7T ) 2 : the decomposition of N2/D4 into partial fractions 

N 2 (r) 71 72 

H hex. (3.37) 



with relations like 



/^(r) r — ct\ r — a 2 



Re £ 7, = 0, 

i=l,2 

2Re ^ «i7i = m 2 , 

i=l,2 



Re ^ ai(ai7i - ""4) = (3.38) 



i=l,2 
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shows that the Laplace inverse j(fj) of rN 2 /D± 

7 (/i) = 2Re ota^ (li > 0) (3.39) 

i=l,2 

as fi — > takes on the value 

7(0) = ml (3.40) 

with a slope of 

7 '(0) = 2771^(1 + m^Re ^ «i) • (3-41) 

i=l,2 

This slope is positive and, due to the small value of m,, it is not very sensitive to small 
variations in Re(«i + 0:2). 

Altogether, we obtain for the isovector electric spectral function 

r /=1 (t') = r cont (t') + Er«(t') (3.42) 

i=l 

where the continuous part is given by 

r co „t(/i 2 ) = ^-i{v-2m n ) (3.43) 
2/i 

with ( |3.39| ) and the amplitudes of the discrete pole structures (|3.S|)-( |3TTl| ) are fixed through 
(PI). 

As before it turns out that the quality of the fit is not very sensitive to variations of 
e 2 within (0.2 < e 2 < 0.4). Keeping e 2 fixed at 0.3 the positions of the roots of D A are 
determined as 

a x = -0.550 ±i 0.215 ; a 2 = -0.3 ± i 0.626 ; [m^ 1 ] (3.44) 
with pole locations and amplitudes given in table II. 



r(o 


i = 1 


i = 2 


i = 3 


% = 4 


Ui [m g ] 


4.39 


0.72 


2.60 


2.92 


di [m 2 g ] 


-7.59 


0.009 


-5.59 


3.49 



TABLE II. Positions and strengths of the discrete structures for the isovector electric spectral 



function (in the soliton restframe) 



The stable feature of the resulting isovector spectral function is that there is no significant 
low-lying discrete strength, while the continuous part rises steeply from its finite limit 

L cont (/i 2 - (2m,) 2 ) - ^ (3.45) 

at threshold to a pronounced maximum near 1.5m e (see fig. 4). Again the discrete structures 
above 2m e mainly reflect our specific choice ( |3.35| ). 
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C. Magnetic form factors G° M (t) and G M (t) 

In the purely pionic Skyrme model, the normalized magnetic isoscalar and isovector form 
factors are given by (t — —Q 2 ) 



G\(t) = - I rf 3 r MQr)^!- (3.47) 



Thus, with the convenient supplementary definition 

Biir 
2 j - • ^v-v ) r 

we have 

G° M (t) = ^ jG%{t), G^t) = 6^<%(t) • (3.48) 

and, correspondingly, for the spectral functions T(t') 

rS#(0 = ^ ' r ^ (t ' } = 6 |^ (t,) (3 - 49) 

for timelike t' = fi 2 . 

The Laplace transformation ( |3.3| ) relates T^(/i 2 ) and Tg(/^ 2 ) by 

/if^(/i 2 ) = r f uT^u 2 ) dv dv' (3.50) 

J HO Jfl 

(with fj, Q = 2m v ) . 

For the continuous part of the isoscalar magnetic spectral function we obtain from (|3.26 ) 
and (ET49p 

r° Mcmt (» 2 ) - ^|l^ 7(/i_3mw) (3 ' 51) 

with 7(/x) given by ( p.23| ). At threshold /x — > 3m,,- T^- rises like 

nicont^ 2 ) - ~ (3m.) 2 ) 2 . (3.52) 

The discrete monopole structure near fi ~ 1.6m g which appeared as a stable feature in 
r^(/i 2 ) is turned into a dipole by the derivative in ( |3.49| ). Therefore T° M contains no discrete 
monopole strength. 

From ( |3.49| ), ( |3.50| ), ( |3.39| ) and ( |3.43| ) we obtain the continuous part of the isovector 
magnetic spectral function 

cont {f) = A Re £ ^ (l + - l)e^- 2 ^) . (3.53) 

" i=l,2 "* 

These continuous parts of the magnetic spectral functions are displayed by the dashed lines 
in figs. 3,4. 
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IV. COUPLING TO VECTOR MESONS 



We have accounted for the influence of vector mesons on the form factors by multiplying 
the purely pionic parts G n (i) with 

A(t) = XD v (t) + 1 - A (4.1) 

where D v {t) is just the vector meson propagator. To illustrate its influence on the spectral 
functions we consider the isovector electric case where the pionic part of the spectral function 
is given by ( |3.42| ) 

r^(t / ) = r CO nt(t , ) + Er w (0 (4.2) 
i=i 

with the continuous part r cont shown in fig. 4 and positions and amplitudes of the discrete 
part listed in table II. 

If we denote the imaginary part of the ^-propagator Im Dy = T Q (t') we have for the 
combined form factor (for A=l) 

1 r T e (f) T n (t") 

G(t) = Dv(t)CT(t) = - 2 J ^£1 _ki dt > dt ». (4.3) 
For the spectral function T(t') of G(t) this implies 

r(f ) = I (r«(tf) vj p^dt" + r*(f) vj ^^'') (4-4) 

where V denotes the principal value. The first integral gives the resulting strength distribu- 
tion of the vector-meson pole. It receives contributions from the continuous and the discrete 
part of which are slowly varying functions of t' within the resonance region. At t'=m v 
we obtain the values 

/■ r cont (r) /- Er(»(t") 

J t"-m 2 v ' J t"-m\ V ; 

so that the total strength of the vector-meson pole is increased by a factor of 2.14. 

The second integral in ( |4.4|) modifies the pionic spectral function. The discrete multipole 
structures rw in T 71 of course remain discrete, but receive additional lower multipoles at the 
same position. (Only monopoles remain monopoles, with modified strength). Of particular 
interest is the modification in the continuous part of the spectral function: The t'-dependence 
of the second integral in ( [4.4|) cuts down the tails in r cont (t / ) for larger values of t' and creates 



a sign change near the vector-meson pole. The result is shown in fig.5: For the demonstrative 
purpose we use for the £>-pole the imaginary part (dotted line in fig.5) 

3 

{ ] (m 2 v - 1') 2 + m 2 v e 2 1 1 

with the experimental full width e=150 MeV (and a smooth cutoff near the two-pion thresh- 
old). The dashed line in fig.5 is the continuous pionic part r cont from ( |3.43| ) and fig. 4. The 
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dash-dotted line shows the continuous spectral function r^T -1 ^') as resulting from ( |4.4|) . 
The increased strength of the resonance due to the integrals ( [4.5| ) is evident. Furthermore 
the figure clearly shows that the interplay between the vector meson and the solitonic pion 
cloud leads to appreciable enhancement of the resonance structure on the low-energy side 
and a corresponding depletion on the high-energy side. This is strongly reminiscent of the 
Frazer-Fulco unitarity constraint which enforces a similar behavior for the isovector spectral 
functions above the two-pion threshold |T3,H0j. Notable from the figure is also the finite 



limit ( ggg ) of T^it^Aml) = hm^/A at threshold. 

The actual value of A used in sect.l for the form factors was A=0.75. Therefore the total 
spectral functions are obtained as the corresponding superpositions of Ar(t') from ( |4.4| ) and 
the purely pionic (1 - A)r^(t') of ( ggg ). This final result for Y K E \t') is shown by the full 
line in fig. 5. 



V. RELATIVISTIC CORRECTIONS 



The Lorentz boosts ( |2.16 ) and ( 2.17| ) which relate the form factors as evaluated in the 



soliton restframe and in the Breit frame can be directly transferred to their respective 
spectral functions. If we denote by T^ nr \t') the spectral function as evaluated in the soliton 
restframe, we have from ( |2.16| ) 



G E (t) = - ^ f 

n Jto v — r/7 



{m V) 



dt' 



f Le W7 j dt + G^ r \-AM 2 ) (5.1) 



with 



Thus, with 



7T Jt t — t 



7 2 = 1 l — for timelike t > . (5.2) 

' AM 2 v ' 



G E (t) = - J^- dt + const. (5.3) 



71 Jto t — t 



T E (i) = Tt\t/f). (5.4) 



we have 



The additional constant which appears in ( |5.1| ) and ( |5.3| ) is the undesirable consequence of 
( gig ) for t -> oo 

G E (t -> ±oo) -> G { E r \-AM 2 ) (5.5) 

which we have met already in ( p.!8|) . 

The integration interval in (|5.3|) is determined by the threshold of T^ E r \ i.e. t/'y 2 = to, 

or 
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The upper limit is the physical threshold for soliton-antisoliton creation which follows from 
(U) for 7 2 -> 0. 



For the magnetic form factors the boost to the Breit frame ( 2.17|) leads to 

i r lM 'r M (t) j: 



G M(t) = L r ^dt (5.7) 

7T Jt t — t 



with 



M«) = ^\t/f) . (5.8) 

In this case there remains no additional constant term in fl5.7|) such that Gjw(t) satisfies 
Gm{1 — * ±oo) — ► 0. The factor I/7 2 in the relativistically corrected magnetic spectral 
functions diverges as t — > AM 2 . This amplifies the oscillations of the non relativistic (t') 
for t' ^> m 2 which are not well determined by the fit (through £2) and thus should not be 
taken too seriously. 



The requirement that the threshold to i n (5.3) should appear at multiples of the physical 
pion mass to = {nm^) 2 (n = 2 or 3 for isovector or isoscalar spectral functions) can only be 
met if in the soliton restframe an effective pion mass /u, n is employed: 

x2 (nmj 2 AM 2 

t = (n^) 2 = A— ±L 5.9 

AM Z — [n m-v) 

This is an extension to timelike to and to of the usual 'reduced mass' formula 

2 [nm^f AM 2 

(n U n ) =——77, 7 (5.10) 

v h ' AM 2 + (n m n ) 2 v ; 

for spacelike momentum transfer. This effect, however, is very small. 

The main effect of the transformations ( |5.4j) and (|5.8|) is that the point t' = 00 is mapped 
onto the threshold t = AM 2 , i.e. the nonrelativistic spectral functions are squeezed into the 
finite interval t < t < AM 2 . For the continuous parts of the pionic spectral functions 
this is shown in the lower sections of figs. 3,4. Of course, the low-energy region is almost 
unaffected, especially the region around the vector meson resonance (after inclusion of the 
£>-pole) remains essentially as shown in fig. 5. But the positions z/j of the higher discrete 

structures listed in tables I, II are shifted to z/j = i^/yl + v 2 /AM 2 . 

If we here use for M the physical nucleon mass M^=940 MeV in order to relate these 
shifted positions to the actual physical N-N— threshold we obtain the values listed in table 
III. Although we cannot really expect any close correspondence to the actually observed 
isoscalar ($(1020), w(1420), w(1600), $(1680)) and isovector (^(1450), g(1700)) resonances 
it is satisfactory to find the calculated structures in this energy range. (The isovector 
structure at 0.72 m g appears with almost vanishing strength 0.009 in table II and may 
indicate that the continuous part expects small modifications near 5 m n .) 
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isoscalar 


isovector 




i = 


i = 1 


i = 2 


i = 3 


i = 1 


i = 2 


i = 3 


i = 4 


Ui [m e ] 


1.61 


5.14 


1.93 


2.09 


4.39 


(0.72) 


2.60 


2.92 


h [MeV] 


1034 


1695 


1165 


1221 


1640 




1369 


1440 



TABLE III. Positions of the calculated discrete structures: Vi, in the soliton restframe (taken 



from tables I, II), and z/j, in the Breit frame (as obtained from z>j = Vij J\ + v 2 /AM 2 with M=Mjy). 



VI. SUMMARY 



In view of the widespread prejudice that baryons considered as solitons in effective 
mesonic theories provide only poor approximations for the e.m. nucleon form factors within 
a very limited range of Q 2 we demonstrate here that even in the most simple soliton model 
the essential features of form factors observed over three orders of magnitude in momentum 
transfer are naturally reproduced. The analysis shows that in order to achieve such a result 
it is necessary to incorporate three basic ingredients: an extended object, partial coupling to 
vector mesons, and relativistic recoil corrections. We have used here the standard skyrmion 
as the extended object, and (for simplicity) the same vector meson propagator in both 
isospin channels. Apart from the Skyrme constant e which is adjusted to the Delta- nucleon 
mass split the model contains only one parameter of direct relevance to the form factors: 
the amount A which measures the admixing of photon- vector meson coupling. We use it to 
bring the e.m. proton radii close to their experimental values. This is sufficient to keep the 
form factors G P E M and G^ very close to the standard dipole up to almost 1 GeV/c. 

Naturally, the typical shape above Q 2 ~ is sensitive to the recoil corrections. We 
employ the boost transformation to the Breit frame suggested by Ji [I1J which pulls the point 
Q 2 = AM 2 (in the restframe) to Q 2 = oo (in the Breit frame). Thus the behavior of the form 
factors for Q 2 — > oo is, of course, very sensitive to the value of M, the kinematical mass of 
the particle represented by the soliton. In tree approximation M is the classical soliton mass, 
but it is known to be subject to large quantum corrections. The boost to the Breit frame has 
the undesired feature that it violates super convergence unless AM 2 coincides with a zero in 
Gm (or a double zero in Ge)- It should be noted, however, that the boost transformation is 
not an exact result, but its derivation relies on neglecting all commutators between center- 



of-mass and total momentum [Oj]. Thus it appears as a reasonable procedure to subtract 
the very small undesired terms from the boosted form factors. We have shown here that up 
to 10 (GeV/c) 2 their effect is qualitatively unimportant and the typical features observed in 
G P M (which is the only form factor known sufficiently well up to high Q 2 ) clearly emerge from 
the model. Due to these difficulties connected with the treatment of recoil corrections the 
predictive power at tree level is poor for the high-Q 2 behavior. This flexibility is unfortunate 
because it prevents definite conclusions and conceals possible evidence for importance of 
QCD effects. From this point of view it would be highly desirable to learn more about loop 
and rotational corrections to the form factors. 

The continuation to timelike t leads to quite stable results for the spectral functions in 
the regime from the 2- or 3-pion threshold to about two rho masses. Beyond this region well- 
known ambiguities prevent definite results. Especially the onset of the continuous part of 
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the spectral functions at threshold can be reliably determined and there are strong analogies 
to the results imposed on dispersion theoretic approaches by the unitarity constraint. On 
the other hand, the region near the physical nucleon-antinucleon threshold t = 4M 2 is most 
poorly defined, because through the boost transformation it is an image of the nonrelativistic 
form factors near Q 2 — > oo. Both, Gffi and G^ r \ vanish in this limit, so the possibility of 
finite results depends on the power of (diverging) factors (1 — t/4M 2 ) -1 introduced through 
the boost transformation. This does not seem a reliable way for investigating this regime; 
it might prove more promising to try to extract the timelike form factors above the N-N 
threshold directly from soliton-antisoliton configurations. 
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FIG. 1. Proton magnetic and electric form factors G V M and G V E (divided by the standard dipole 
Go) as obtained from the Skyrme model with e=4.25, (soliton mass M=1648 MeV), and vector 
meson coupling A=0.75 in ( |2.13 ). The dashed lines show the same form factors after subtraction of 
the square brackets in (^Tl8|) , ( |27T9|) . The data are from the compilations of [jl4|] (open circles), |3^] 
(open squares), (full circles), and [^8| (full triangles). 



FIG. 2. Neutron magnetic form factors G% (divided by the standard dipole Gd) and neutron 
electric form factor G E (full lines) calculated as in fig.l. The data for are from [^f| (full 
circles), ]27| (diamonds), Q| (triangles), and |28| (square). The open circles represent older data 
from |29|-^4|. The data points for G E are from |35| for the Paris potential (full circles) and 
from [36] for Feshbach-Lomon wave functions (open circles). The dash-dotted line for G E is the 
Galster fit [36], the dotted lines are fits from f3f| deduced for three different potentials. 



FIG. 3. The continuous parts r cont of the isoscalar electric (full line) and magnetic (dashed 
line) spectral functions for the skyrmion with e=4.25. In the upper part they are shown in the 
soliton rest frame, in the lower part they are boosted to the Breit frame (M=1648 MeV). 



FIG. 4. Same as fig. 3 for the isovector spectral functions. 



FIG. 5. The continuous part of the isovector electric spectral function after inclusion of the 
£>-resonance. The dotted line is the imaginary part of the £>-propagator (with width 150 MeV); 
the dashed line shows the continuous pionic part r cont from ( 3.43 ) and fig.4; the dash-dotted line 
shows r^T -1 ^') as resulting from (|4.4|) (i.e. for A=l). The full line is the total result for A=0.75. 
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